Regression modelling in hospital epidemiology: a statistical note

Barnett and Graves [1], in their commentary on our report recently published in Critical Care [2], suggested that time-discrete methods should be used to address time-dependent risk factors and competing risks. In this letter we comment on two statements by those authors. First, Barnett and Graves claim that, 'An alternative method to the competing risks model is a multistate model.' In fact, a multistate model is not an alternative to modelling competing risks, but a competing risks model is an example of a multistate model. This is explained in the tutorial by Putter and coworkers [3]. However, competing risks only model the time to first event and the event type (for example, time to noso-comial infection [NI]) or discharge/death, whatever comes first. To model subsequent events also, more complex multi-state models are needed. Barnett and Graves give an example in which discharge/death events after NI are also modelled. However, such a complex multistate model is implicitly used in a competing risks analysis when time-dependent risk factors are included. For example, in our report we also analyzed discharge and death as competing events and NI as a time-dependent covariate. Such a model analyses discharge/death events after NI. Second, Barnett and Graves claim that pooled logistic regression has some advantages over Cox regression. In our intensive care unit setting, if data are collected day by day then both models are technically identical survival models and provide identical results, even though the Cox model is less restrictive. The Cox model is a very flexible regression model with potential extensions. Usually, hazard ratios are assumed to be constant, but the assumption of proportional hazards can be relaxed such that the regression coefficients may vary with time (see Martinussen and Scheike [4]). Time since NI can be included as an additional covariate in this model. Random effects may also be studied (see the discussion of frailty models in the book by Hougaard [5]). However, in this context it should be noted that standard logistic regression (one patient-one record) does not model time to event and might not appropriately include time-dependent risk factors. In accordance with our primary report [2], Barnett and Graves [1] highlighted the need to take into account both time-dependent risk factors and competing risks when the impact of nosocomial infections on intensive care unit death/ discharge is the focus of interest. In addition, they highlighted important issues that would be worth studying …